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ABSTRACT 

Faraday rotation (RM) probes of magnetic fields in the universe are sensitive to cosmological and 
evolutionary effects as z increases beyond '--^l because of the scalings of electron density and magnetic 
fields, and the growth in the number of expected intersections with galaxy-scale intervenors, dN/dz. In 
this new global analysis of an unprecedented large sample of RM's of high latitude quasars extending 
out to z ~3.7 we find that the distribution of RM broadens with redshift in the 20 — 80 rad m~^ range 
range, despite the (1 +z)~'^ wavelength dilution expected in the observed Faraday rotation. Our results 
indicate that the Universe becomes increasingly "Faraday-opaque" to sources beyond z ~ 2, that is, 
as z increases progressively fewer sources are found with a "small" RM in the observer's frame. This 
is in contrast to sources at z <1. They suggest that the environments of galaxies were significantly 
magnetized at high redshifts, with magnetic field strengths that were at least as strong within a few Gyr 
of the Big Bang as at the current epoch. We separately investigate a simple unevolving toy model in 
which the RM is produced by Mgll absorber systems, and find that it can approximately reproduce the 
observed trend with redshift. An additional possibility is that the intrinsic RM associated with the radio 
sources was much higher in the past, and we show that this is not a trivial consequence of the higher 
radio luminosities of the high redshift sources. 



Subject headings: galaxies: high redshift 
data analysis 

1. INTRODUCTION 
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The strengths of interstellar and intergalactic magnetic 
fields at earlier epochs have important implications for 
galaxy and structure evolution (e.g. Mestel & Paris 1984; 
Rees 1987), the propagation of ultra-high energy cosmic 
rays (Sigl, Miniati, & EnBlin 2003, 2004; Armengaud, Sigl 
& Miniati, 2005; Dolag et al. 2005), and the feedback 
of magnetic energy into the intergalactic medium by stel- 
lar winds and early massive black holes (Kronberg et al. 
2001). 

Faraday rotation of distant polarized radio sources is 
one of the few available measurables to detect and probe 
extragalactic magnetic fields. For a cosmologically distant 
polarized source at redshift Zs it is defined, in units of 
rad/m^, as: 



RM{zs) = 



= 8.1 ■ 10" 



ne{z)B\\{z) dl 

(1 + 2)2 ^ 



dz. (1) 



includes any metagalactic and/or local universe RM con- 
tributions that might exist on large angular scales; (2) A 
component arising from intervening discrete clouds, e.g. 
galaxies, and/or a diffuse medium along the line of sight. 
The latter includes filaments of cosmological Large Scale 
Structure (LSS). The invervening galaxy system compo- 
nent should depend, statistically, only on the intergalactic 
path length traversed and not on direction. The third 
(3) is an "intrinsic" component from magnetised plasma 
associated with the distant radio source and its immedi- 
ate environs, (RRMintr), which may depend on source- 
intrinsic properties and which may also evolve cosmolog- 
ically; Finally (4) there are measurement errors, which 
ideally should not depend on Z, &, or z. 

Detailed RM images of individual quasars between z ~ 1 
and 2 ~ 2 with companion absorption line data have es- 
tablished clear examples of case (2) e.g. PKS 1229-021 
(Kronberg, Perry & Zukowski 1992), and of case (3), e.g. 
3C191 at Zs — 1.95, which contains intrinsic RM varia- 
of order 2000 rad m"^ (Kronberg, Perry 



The RM describes the change in polarization angle Axo 
with respect to a change in wavelength squared AAq due 
to the presence of a magnetized medium (the subscript 
indicates observer's frame). In Eq. (1) the free electron 
number density, rig, is in cm~^, _B||, in Gauss, is the line 
of sight component of the magnetic field and dl/dz, in 
parsecs, is the comoving path increment per unit redshift. 
In general the total RM of a given radio source will be 
a sum of several different components: (1) a "smooth" 
Galactic component, defined as SRM, that may be as- 
sumed to vary with I and b on angular scales that are 
larger than the typical inter-source separation. This also 
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tions at ^ 

& Zukowski 1990). These studies probed the intervening 
Faraday-active gas by analysing both optical absorption 
lines and Faraday rotation images. More recently, two- 
dimensional RM images from larger samples of resolved 
high-z quasar radio maps (Athreya et al. 1998, and Carilli 
et al. 1997), also show similarly high Faraday rotations at 
2 < z < 4, which these authors interpreted as indicating 
that the intrinsic RM dominates, i.e. case (3). 

Redshift-dependence of the SRM-corrected RRM(z) 
(RRM = RM - SRM) could be due to one or both of (2) 
or (3) above, and in each case is subject to the (1 + z)^^ 
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watering-down effect in equation (1). Earlier attempts 
were made to detect a z-dependence of quasar RM's (Rees 
& Reinhardt 1972, Kronberg & Simard-Normandin 1976) 
using RM data on samples of extragalactic radio sources. 
All of these showed some evidence for an increase in the 
observed RM of quasars at z > 1. With somewhat bet- 
ter RM data, and with optical absorption line data for 
some quasars, it was found that high column density op- 
tical and HI absorption at intervening redshifts correlated 
with higher levels of observed RM (e.g. Kronberg & Perry 
1982, Kronberg, Perry & Zukowski 1992, Oren & Wolfe 
1995). This allowed some first estimates of magnetic field 
strengths in distant intervening galaxy systems. 

Welter et al.(1984) found a clear growth of the overall 
width (variance) of |RRM|(z) in a 116-quasar RM sam- 
ple, recently confirmed in a smaller sample by You et al. 
(2003). Welter et al. also developed mathematical frame- 
works for connecting RM correlations to components (2) 
and (3). They tentatively favored (2) over (3) at redshifts 
up to z ^ 2. A framework has also been developed by Ko- 
latt (1998) for an RM-based probe of the primordial mag- 
netic field spectrum, related to cases (2) and (3). More 
recently, optical galaxy and RM data were combined to 
undertake the first magnetic field probe of LSS filaments 
in the local universe (Xu, et al. 2006). Magnetic fields in 
the cosmic voids of LSS have not yet been detected. 

This paper analyses an optimized subset of a new, much 
expanded sample of 900 extragalactic RRM's with mea- 
sured redshifts up to z ~ 3.7. These also have improved 
determinations of the (smoothed) Galactic and local uni- 
verse foreground SRM contribution. We focus on the over- 
all distribution of RRM's as a function of redshift as a 
diagnostic of magnetic fields in high redshift systems, in 
particular the little-explored z - dependence of |RRM| at 
|RRM| < 100 rad m~^. This is distinct from previous in- 
vestigations, summarized above, of the overall width, or 
variance of the RRM distributions. 

The rest of the paper is organized as follows. In Section 
[2 the RM dataset is briefiy described and in Section [3] we 
discuss the selection and optimization of the RM sample, 
as well as tests for unwanted selection effects. The key 
evidence for a redshift dependence of the RRM distribu- 
tion is presented in Section [H We interpret the observed 
effects in Section [5j In general terms, an increase in the 
width of the RM distribution requires the presence of sub- 
stantial amounts of magnetized plasma at high redshift. 
As one hypothesis, we construct a model in which the RM 
arises due to the presence of magnetized clouds traced by 
Mgll absorption systems (|5.ip . In section 6 we discuss 
the complexities of detecting an all-pervading, intergalac- 
tic magnetic field in the presence of other extragalactic RM 
contributions. Our conclusions are summarized in Section 
7. Throughout the paper we assume a concordance cos- 
mology with Ho = 70 km s^^Mpc^^ = 0.25, f^A = 
0.75. 

2. THE DATASET 

The expanded RM data include new linear polarization 
measurements by one of us (PPK), which favored quasars 
at larger redshifts, plus additional data from the litera- 
ture. Our new sample of 901 quasars and radio galaxies is 
sufficiently large that we are able to choose an optimized 



subset (by Galactic l,b location) of 268 objects distributed 
up to z ~ 3.7. This subset exceeds the previous all-sky 
116 quasar sample of Welter et al. (1984) by a factor of 
2.3. In addition, the average quality of RM measurement 
is improved, as is also the Galactic foreground correction, 
which used a new expanded sample of 1566 extragalactic 
radio source RM's, mostly at |6| > 5°. 

Accurate subtraction of the foreground Galactic SRM 
is important for obtaining the purest possible isolation 
of the extragalactic RRM. The SRM was estimated using 
a Bayesian, Gaussian process formulation applied to the 
larger all-sky sample of 1566 extra- galactic radio source 
RMs (Short, Higdon and Kronberg, 2007a, b) which gives 
an SRM estimate for any {I, b) location. 

In that paper, a fairly general Gaussian process (GP) 
model for the surface of the unit sphere S with great circle 
distance metric d{si, S2) is derived by taking a collection 
of uniform, regularly spaced knot locations Wi,...,wj, 
and assigning to each of these locations a knot value 
xi,...,xj which are assumed to have iid A''(0, [JAa,]"-^) 
distributions, where iid denotes "independent and identi- 
cally distributed." Convolving these knot values with a 
simple smoothing kernel fc(-) then results in the GP model 
J 

u{s) = ''^Xjk(d{s,Wj)), seS. (2) 

j=i 

Figure [T] shows an example where S is the unit circle. As 
J — !■ 00, the process u(s) quickly converges to a stationary 
GP. 

The convolution kernel k(-) is taken to be a normal den- 
sity whose width is to be estimated. A recursive tessella- 
tion algorithm is used to distribute J = 2562 knots over 
the unit sphere, giving a neighboring knot-to-knot distance 
of approximately 27r/80. The spatial SRM field u{s) is 
given by ([2]) which requires that the knot values xj and 
the kernel width be estimated from the radio source RMs. 
The N — 1566 observations taken at locations Si, . . . ,sn 
are modeled as 

Y{s,) = u{si) + e,, i^l,...,N 

where the errors are independent with A^(0, [wiAe]"^) dis- 
tributions. The Wi's, which modify the error precision, 
account for the possibility that certain observations have 
been altered by individual source RM anomalies, inter- 
veners, very small scale local effects of the Milky Way, 
etc. When oji < 1, the observation has been altered and 
is downweighted for the purpose of estimating the SRM 
field; when LUi = 1, the observation is assumed to be free 
of any of these altering effects. 

From the analysis of Short et al. (2007a), the propor- 
tion of the observations identified as altered {uji < 1) is 
23%. The resulting posterior distribution from this model 
formulation, which accounts for (a) uncertainty regarding 
the spatial dependence in the SRM field, (b) the classifica- 
tion of observations as altered or not, and (c) observation 
noise, was sampled using Markov chain Monte Carlo. The 
model formulation is described in detail by Short, Higdon 
& Kronberg (2007a). Compared with previous methods 
that averaged the neighboring unedited RM(Z,6), and then 
iteratively remove "outliers" (e.g. Simard-Normandin & 
Kronberg 1980 and other papers since then), the more 
statistically formal method used here is free of ad hoc cri- 
teria, for example on the decision of when to delete an 
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Fig. 1. — A realization of a process convolution model on the unit circle. The left panel shows knot locations and values (positive and 
negative), along with a smoothing kernel. The right panel shows the resulting Gaussian process realization obtained by a convolution of the 
knot values with the smoothing kernel. 



outlier. Another important advance is that our model is 
with the data to produce the best estimate of the kernel 
width, rather than having it estimated, or guessed a priori. 

Because this paper focuses on links between source red- 
shift and RRM, it is important to search for possible bias 
effects of Galactic sky location on the z-distribution of 
the sources. As expected, the Galactic {l,b) locations of 
the RRM's do not correlate with source redshift. Sim- 
ilarly, neither the amount of foreground RM correction, 
SRM, nor the uncertainty in the SRM show any trend 
with source redshift, after we isolated the optimum (Z,6) 
zones as described below. 

3. OPTIMIZATION OF THE DATASET AND INSPECTION OF 
BIAS EFFECTS 

To ensure that other biases have not been introduced 
in the data selection we address three important issues 
that could potentially affect our statistical analysis. First, 
although the subtraction of the Galactic foreground is op- 
timized, the quality of the SRM subtraction removal is 
{l,b) dependent, being generally less accurate at the lower 
Galactic latitudes. This leads to Galactic zone— based ad- 
mission criteria for our sample, elaborated in section [3.11 

Secondly, our sample of radio sources is, typically, 
flux limited so that higher redshift objects tend to have 
have higher luminosities. Thus we also investigate the 
RRM— Liadio correlation to test for possible luminosity se- 
lection effects that could prejudice our search for a RRM 
- redshift relation. 

Thirdly, errors in the individual RM measurements (sep- 
arate from the Galactic foreground correction uncertainty) 
could conceivably exhibit a systematic redshift depen- 
dence, e.g. if the polarized signals were much weaker at 
larger redshifts. On investigation we find no significant 
correlation with redshift of the uncertainty of the individ- 
ual RM determinations. 

3.1. Optimal isolation of true extraglactic RM's 

To further minimize the uncertainties due to the sub- 
traction of the Galactic component, SRM, we examined 
how the form of the overall RRM distribution varies as a 
function of the {l,b) region on the sky. We find that the 
width of the distribution N(RRM) for the smallest RRM's 
decreases steadily as we raise the lower limiting latitude 
from |6| = 10°(lower boundary) up to |6| w60°, but there- 
after asymptotes to a stable value. In two dimensions, 
we find that the boundary can be lowered at some longi- 
tude zones without increasing the Galactic dispersion. An 
optimal subsample of 268 objects for which the influence 



of SRM is minimized was thus obtained by accepting all 
sources at |6| >60° as above, plus those having b > 45° and 
; = 150 - 360°, and those with b < -50° at ^ = 120 - 180°. 

3.2. Checks for luminosity and other unwanted 
systematic effects on RRM(z) 

Given that the radio luminosities of the sample span 
about six orders of magnitude, we can test whether the 
changing mean radio luminosity with redshift is likely to 
be significant. For the 268 radio source sub-sample we 
compared the rest-frame ( "/c-corrected" ) 2.7 GHz radio 
luminosity, L2.7, with the maximum rest-frame RRMjntr- 
The latter assumes that all the RRM comes from mag- 
netized plasma in the vicinity of the source, that is we 
multiply the observed RRM by (1 + Zs)^ (eq. [1]). 

At z < 0.2, where neither cosmological evolution nor in- 
tersection of the line of sight with an intervenor are likely, 
the radio luminosities cover more than five orders of mag- 
nitude (Figure However, we find no evidence for any 
correlation between |RRMinti | and L2.7 (see Figure^]). To 
quantify this statement we parameterize the distribution 
in RRMjntr using a Lorentzian function 

/(RRM;F/2)^l ^^^^^,^/;^^„ (3) 

where r/2 is the width at half maximum, and we test for 
a r/2 - L2.7 relation with 

r/2 = u)o(L2.7/102^WHz-i)T. (4) 

The best fit parameters for this relation at z < 0.2 are: 
wq — 36^24'"ad/m^ and 7 = — O.IIIq :[q, confirming that 
there is no effect of L2.7 on |RRMi„tr|. The right hand 
panel of Figure O shows this graphically. Sources with 
lines of sight used to test for a correlation are represented 
as black dots. The blue stars and red squares are for lines 
of sight from sources at higher redshift. Due to the (l-t-z^)^ 
transformation they produce much larger |RRM| values. 

In conclusion, our test confirms the lack of any obvi- 
ous L2.7 - RRM relation. Note, however, that we cannot 
exclude the possibility that the highest redshift and most 
luminous sources could be affected by strong evolution- 
ary effects due to the presence of enhanced magnetic fields 
when compared to their low redshift/low luminosity coun- 
terparts. 

4. ANALYSIS OF THE RRM-REDSHIFT BEHAVIOR 

Because of the statistical nature of the RM measure- 
ments and the present limited availability of supplemen- 
tary optical data on individual lines of sight, the most 
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Fig. 2. — Left panel: Radio power at 2.7 GHz against redshift for the sample used in our analysis. The solid line marks the subsample 
(z < 0.2) used to test for a correlation between luminosity and intrinsic RRM. Right panel: Luminosity against RRM values transformed 
to the radio source rest frame (multiplied with (1 + Zs)^)- The empty circles are for the subsample at z < 0.2 that was used to test for a 
correlation between luminosity and RRM, the stars are those with 0.2 < z < 1.0 and the solid squares are in the range 1.0< z < 3.7. 



powerful diagnostic that we have is the /orm of the ob- 
served distribution N(RRM, z), as a function of redshift. 

To illustrate for the case of a source-intrinsic RM, if this 
were 2;-independent the observed RRM should monotoni- 
cally decrease with increasing source redshift, due to the 
strong term (eq.([T|)). This means that any varia- 

tions in N(RRM, z) at large z will arise from a competition 
between genuine evolutionary effects and the (1 + z)^^ re- 
duction. Thus, an interplay of different effects at 2 > 1 
might be expected. 

On the other hand, in the case of intervening systems 
at low redshift their contributions to |RRM| will be sta- 
tistically invisible until there is a significant probability of 
intersecting an intervenor, at which point they will rise. 
Then at some sufficiently large 2, the effect of incremen- 
tal RM intersection depth may be overwhelmed by the 
(1 -|- z)"^ term, which would tend to flatten any increase 
in the observed RRM. It may subsequently increase again 
if a sufficiently strong evolutionary increase sets in again 
at still larger redshifts. These possibilities set the context 
for our analyses in Section [5] 

In the following we use a non-parametric approach, 
counting fractions of lines of sight ji below a varying 
threshold RRMi as a function of z (see Welter et. al 1984). 
We then apply a Kolmogorov-Smirnov test to examine the 
statistical significance of the redshift dependence in our 
data. 

4.1. Evidence for Redshift Evolution 

Figure [3] shows the observed distributions of RRM for 9 
different redshift bins, each containing about 29 sources. 
At redshifts below z 1 the distributions are character- 
ized by a sharply defined mode at RRM w plus a broader 
RRM component. Beyond this redshift, and especially at 
z > 2, the low-RRM component tends to get redistributed 
to larger RRM values. 

This effect is shown quantitatively in Figure [H where 
the fractions of lines of sight fi below a varying threshold 
RRMi are plotted as a function of redshift. The quan- 
tity /20J representing the fraction of |RRM|'s below 20 
rad m~^, decreases significantly with redshift, from 72 % 
in the lowest two redshift bins {z = 0.08) to 39% in the 
highest redshift bin (z = 2.34). Similar, but progressively 
weaker trends occur for /40 and /loo respectively. The 
highest redshift bin shows the smallest value of fi for all 
three thresholds. These trends indicate a clear evolution- 



ary pattern which is most apparent at |RRM| < 40 rad 
m~^ in the observer's frame. 

We confirm the reality of this redshift dependence with 
a Kolmogorov-Smirnoff test. The entire sample was di- 
vided into two, below and above Zb and we compare the 
two normalized cumulative distributions of the absolute 
value of the RRM, N(|RRM|) in Figure [5] (left panel). For 
Zb ^ 1.8, the samples are different at the 99% significance 
level. The RRM at which the maximum perturbation in 
the N(|RRM|) occurs is typically between 10-25 rad m^^ 
(Figureini right panel). This shows that the clearest evolu- 
tionary signal comes from a broadening of the low- 1 RRM | 
peak, rather than from outliers in the high RRM tails of 
the distribution. 

The "migration" of some sources near |RRM| = to 
wider wings works oppositely to the expected (1 -I- z)^^ 
decrease. It clearly demonstrates that there is a global 
evolutionary effect in Faraday rotation. More precisely, as 
z increases up to and beyond 2, there is a gradual broad- 
ening of the distribution of RM's and a corresponding "de- 
population" of the lowest |RRM| bins, especially those < 
|20| rad m~^. In effect, at higher redshifts the Universe 
becomes progressively "Faraday opaque" as fewer sight- 
lines are able to "escape" an enhanced Faraday rotation 
at the higher z's. The extra Faraday rotation is produced 
at z > 1 and is consequently greater in the Faraday rotat- 
ing rest frame than what we observe, by a factor of (1-t-z)^. 
For example, |RRM|'s of 20 to 100 in Fig. [3] become 80 
to 400 rad m'^ at z = 1 and 360 to 1200 rad m^^ ^j^g 
Faraday rotation originates at z = 2.5. 

Another striking illustration of this redshift effect is 
shown in Figure [SI which compares the normalized cumu- 
lative counts N(< (I RRM I) vs. |RRM| above and below 
Zb = 1. Evolution in the high- |RRM| tails remains to 
be better specified in future, larger samples at the high- 
est redshifts, and we do not attempt to quantify it here. 
What is clear is that a significant evolution in the observed 
|RRM| at modest RRM's begins to set in beyond z ^ 1.0. 

5. QUANTITATIVE MODELLING TO ESTIMATE MAGNETIC 
STRENGTHS IN HIGH REDSHIFT SYSTEMS 

In the following two subsections we take different ap- 
proaches to the analysis and interpretation of the data. 
Each draws from the analysis framework developed in Wel- 
ter et al. (1984). 

First, in Section 15.11 we relate these results to recent 
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Fig. 3. — |RRM| histograms of 268 lines of sight for different redshift bins having approximately equal numbers of lines of sight per bin. 
Redshift increases from the upper left to the lower right panel. The distribution in the highest redshift bin (z > 1.8) is characterized by a 
significantly decreased "peak" near |RRM| = 0, and a corresponding significant broadening of the "low-|RRM| " peak to higher |RRM's|. 
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Fig. 4. — Fraction of lines of sight fi below a threshold RRMi as a function of redshift. Circles, triangles and squares show f20, f40 and fioo 
respectively. Errors are calculated by randomly drawing lines of sight and calculating the r.m.s. of the resulting fi distributions. 
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Fig. 5. — Left panel: Significance level (SL) of the KS test that the distributions of the RRM's below and above zt are not drawn from 
the same parent distribution. The black solid line is with no |RRM| cut and the blue line is for lines of sight with |RRM| < 200 rad/m^. 
The dashed dotted lines give the fraction of lines of sight which are below zi, (right-hand axis). Right panel; The RRM value at which the 
normalized cumulative distributions, N{RRM) from the KS-test most differ, as a function of Zj. 
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Fig. 6. — A comparison of the normalized cumulative counts of N(< |RRM|) vs. |RRM| shown separately for sources having z < 1 (upper 
curve) and z > 1 (lower curve). The clear separation between these two curves shows that excess RRM's introduced in the higher redshift 
subset are typically in the range |RRM| ~ 20 to ~ 80 rad m"'^ in the observer's frame. 



Table 1 

Parameter values for different RRM models. 



model 


c^noise (rad/m^) 


r,;„tr/2(rad/m2) 


o-cioud (rad/m^) 


no intervening systems 








Mgll systems with Wr > 0.02A 




7+6 
' -4 


60+'fi 


Mgll systems with Wr > 0.3A 


9ti 


7+6 
' -4 


1 1 ^+45 



quasar optical absorption line data, given earlier evidence 
for an RM-absorption line association (e.g. Kronberg & 
Perry 1982, Welter et al. 1984, Oren & Wolfe 1995) and 
the a priori expectation that high column density absorp- 
tion line systems will have an effect on RRM. We make 
model predictions of the N(RRM, z) behavior based on 
previous Mgll absorption line studies of quasars up to 
z ~2.3. In order to keep the number of free model param- 
eters appropriate to the current state of the RRM data we 
assume no local cosmological evolution in the RRM inter- 
venor systems, and we also exclude the very high RRM 
outliers from the analysis. We will show that, even in the 
absence of z-evolution and other contributions, Mgll in- 
tervenor systems can have a recognizable influence on the 
observed RRM behavior described in §4. 

In Section 15. 2[ we interpret the broadening in the dis- 
tribution of N(RRM) at z > 1.5 to draw general conclu- 
sions about the strength of early universe magnetic fields 
in galaxy systems up to z 3.5. 



5.1. RRM Intervenor model based on Mgll absorber 
statistics 

In this section we explore the possibility that Mgll ab- 
sorption systems are magnetized: Given their redshift dis- 
tribution from QSO absorption line studies we investigate 
whether they can explain the observed statistical proper- 
ties of the RM data and if so, what the implications are 
for their magnetic properties. Due to the poor statistics in 
the high-RM-tail of the distribution here we consider only 
lines of sight with an observed RRM value smaller than 
|200| rad m^^. Following Welter et al. (1984) we calcu- 
late the probability distribution function P(RRM,2s) for 
an observed RRM value of a source located at redshift Zg, 
as 



Umax 

P(RRM, z,) = qn(zs)Pn(RRM, z,). (5) 

Here P„ is the normalized probability distribution func- 
tion of RRM for a line of sight to a source at redshift Zg 
passing through n intervenors. In addition, qn{zs) is the 
probability of having n such intervenors along the line of 
sight which is given by Poisson's statistics: 

9„(z,) - (n!)-ij.re-''% (6) 

with 

i^s = / -r- dz (7) 

Jo dz 

the mean number of intervening system out to z^ cal- 
culated using the Mgll absorber distribution, dN/dz, de- 
scribed below. In addition to the effects of intervenors 
we also allow for contributions to the observed RRM from 
an intrinsic component and measurement errors, the latter 
dominated by uncertainties in the removal of the Galac- 
tic contribution. As a result, P„(RRM,Zs) is given by the 
convolution of the probability distribution functions asso- 
ciated with each individual component, namely 

P„(RRM,Zs) = Pnoise * ^ n,interv i^s) * Pintri'^s): (8) 

where for both the intervenor and intrinsic component we 
have explicitly indicated the redshift dependence. 

In order to model P„^i„tert)(RRM, Zg) we assume that 
each intervenor can be described as a cloud characterized 
by a number density of free electrons Tie, a size L and 
a randomly oriented magnetic field B with a coherence 
length Ic- Then for each intervenor we can define a prob- 
ability distribution function, Pc/o«d (RRM, z), given by a 
Gaussian distribution with a{z) = acioud{^ + where 
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f^cioud ocUeBlc {L/lcy/^, i.e. has units of RM. The con- 
tribution from n intervening systems is then given by the 
expression 



115 



+45 



Pr, 



,{RRM, zs 



ACn 



dN 

Pcloud{RRM,z) — {z) 

az 



(9) 



which is the convolution (C,i[-]) of n identical clouds dis- 
tributed up to redshift Zs according to dN/dz, and A is a 
normalization factor. Note that expression [S] is valid even 
if Ic ~ L, provided that the magnetic field is randomly 
oriented from cloud to cloud and that the analysis is ap- 
plied to a large number of RRM's associated with different 
lines of sight. 

As for the intrinsic component, typically characterized 
by a non-Gaussian tail, we find it appropriate to use a 
Lorentzian distribution (eq. [3]) which takes better account 
of the outlying RM values. Assuming for simplicity no 
source evolution, the latter can be fully characterized by a 
constant rest frame half width at half maximum, Tintr/'^, 
which translates into an observed T{z) = Tintr (1 + -z)"^- 
Finally, Pnoise, the combination of the observational error 
in RM and the uncertainty in the SRM can be represented 
by a Gaussian width, (Tnoise, which is independent of red- 
shift. 

In the following we carry out two separate analyses, one 
restricted to strong absorbers, i.e. those with an equivalent 
width Wr > 0.3 A, and a second including weak absorbers, 
i.e. those with an equivalent width 0.02 A < Wr < 0.3 A. 
For simplicity all the absorbers are characterized by the 
same adoud, independent of redshift and equivalent width 
(or underlying column density). 

For the weak Mgll absorbers we use the function dN/dz 
obtained by Churchill et al. (1999) who investigated 
HIRES/Keck spectra of 26 QSOs in the redshift range 0.4 
< z < 1.4 . They obtained the result 

^^^^ = (0.8±0.4)(1-Hz)i-™. (10) 

Similarly, for the strong Mgll absoption systems we uti- 
lize the results of Nestor et al. 
0.4 < z < 2.3, find 



(2005) who, for the range 



dz 



= 1.001 (1 
-0.3 



\0.226 



-6.0 



^"pU(i+-)'^J"'"H«(i+^)^' ' ^''^ 

where a = 0.443 and /3 = 0.634. Note that the different 
distributions of weak and strong intervenors, which enter 
the Poisson statistics (O through the parameter i/g defined 
in eq. ([7]), imply different values for the redshift at which 
an absorber is first encountered. 

Using the probability distribution function in eq. O with 
the specific functional form detailed above, we now present 
the results of a maximum likelihood analysis which yields 
the parameters that best describe the data. 

We begin by considering the infiuence of an RRM com- 
ponent due to strong Mgll absorbers, in addition to an 
intrinsic component and Galactic foreground correction er- 
rors. The value of the parameters that produce the best 
fit to the data, in units rad m~^, are 



rmtr/2 

and 
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-1-6 



a - q+* 

'-'noise — '^—2' 

We note that, intriguingly, the typical derived value for 
o'cioud is not very different from an RM observed for a 
typical line of sight through a (low z) spiral galaxy. 

In order to illustrate how our model for P(RRM,Zs) 
describes the data we define the quantile RRMx via the 
relation 

<.+RRMx(z) 

X= P(RRM, z)dRRM, (12) 

J-RRMx(z) 

that is the value such that a fraction X of all the RRM's 
in a distribution P fulfills \RRM\ < RRMx. In Fig. [7]we 
plot both the observed data (points) and the model (solid- 
line) as a function of redshift for the quantiles RRMo.es 
(left) and RRMq.qo (right). 

It is apparent that our simple model with magnetized 
strong Mgll absorbers is capable of reproducing the data 
distribution up to z 2. The redshift dependence of the 
data in this model can be interpreted as follows. At low 
redshifts very few lines of sight pass through an intervenor 
so that P(RRM, Zg) is dominated by the intrinsic and er- 
ror components and varies only slowly with z. As z ap- 
proaches ~ 1, however, the probability of intersecting an 
absorbing system becomes significantly higher. Here the 
role of intervenors starts to kick in and RRMo.es grows 
gently. As we move to sources at progressively higher red- 
shifts, the number of intervenors increases. However, their 
contribution is suppressed by the (1 -I- z)"^ dilution factor. 
This leads to a flattening in the distribution of RRMo.es at 
z ~ 2. Similar reasoning applies to RRMo.go except that 
the few intervenors encountered at low redshifts are able 
to affect this statistical quantity much earlier. 

To illustrate the robustness of our estimate of Udoud, we 
show in Figure [8] contours of constant likelihood as a func- 
tion of the parameters Udoud and Tintr/'^ for two different 
choices of t noise, which will bracket a good fraction of the 
values this parameter can assume. It is apparent that a 
change in Unoise affects Tintr/'^, and that these two quan- 
tities are anticorrelated. This makes sense in that r2nir/2 
and Gnoise dominate at low redshift where, to some extent, 
their role can be interchanged. Importantly, however, the 
uncertainties in the error estimate do not seem to have an 
impact on adoud, which remains at about 3 a above the 
null value. 

When we repeat the same analysis including the con- 
tribution of weak absorbers, the quality of the model pre- 
diction, represented by the dotted line in Fig. [51 worsens, 
particularly for the quantity RRMo.es- This is because the 
number of intervenors encountered below redshift unity is 
substantially higher in the model, leading to a higher ex- 
pected value for RRMo.es- This result may indicate a limi- 
tation in our assumption of a column density independent 
o' cloud for all absorbers. We expect that in attempting 
to reconcile the data with the model, improving on this 
limitation would provide a better solution than postulat- 
ing drastically different magnetic properties for weak and 
strong Mgll absorption systems. Also, establishing a 
= control sample, not available in this investigation, will 
improve future parameter specifications for Mgll interven- 
ers. 
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Fig. 7. — Quantiles RRMq.ss (left panel) and RRMq. 90 (right panel) as a function of redshift are shown for the observed RRM data and 
different models. The solid lines show a model with strong (Wr > O.SA ) Mgll absorbers. The dotted lines shows a intervenor model where 
the weak ( O.O2A < Wr < 0.3A) absorbers are included for comparison. The dashed lines show a model incorporating SRM removal error and 
a non-evolving intrinsic contribution. Each bin contains about 51 lines of sight, and we show the median redshifts. Errors in the quantiles 
are calculated with the bootstrap method, with the Icr confidence interval shown. The derived parameter values are summarized in Table 1. 



Finally, we show a model in Figure [7| (dashed line) in 
which we set adoud = 0, i.e., which has only an unchang- 
ing intrinsic component and a measurement error. The 
dashed lines in both panels of FiglT] show the quantiles 
RRM0.68 and RRMo.90 predicted by this model, with best 
fit parameters, in units of rad m~^, anoise = ISlg and 
rmtr-/2 = 2llg. It clearly fails to reproduce the behavior 
of the observed RRM statistics with redshift. The rea- 
son this model fails so dramatically is that, since the noise 
component is independent of redshift and the RM of the 
non-evolving intrinsic component declines as it 
can only predict a decrease of both quantiles as a function 
of redshift, instead of the observed increase. 



5.2. Implications for magnetic fields in early systems 

Whether due to interveners or other evolution of galactic 
or pre-galactic systems, the rise in the fraction of lines of 
sight affected by Faraday rotation at progressively higher 
redshifts indicates the presence of magnetized gas in such 
systems at earlier epochs. While the intervenor model 
would explain the behavior of the data with the higher 
incidence of magnetized clouds, an alternative possibil- 
ity that we cannot yet exclude is that the RRM in the 
high redshift sources is dominated by an intrinsic com- 
ponent that increases into the past. Note that in this 
case, however, we require substantial negative evolution 
in the Faraday-active medium, i.e. with decreasing 2, in 
order not to over-predict the expected dispersion in the 
observed RRM values at low redshift (see Figure [7]). In 
terms of our model parameters, for example, Ti„tr would 
have to decrease with time. This would follow naturally if 
the Faraday active medium were a magnetized cloud ex- 
panding as a result of being over— pressured with respect 
to the ambient medium. 

In any case we can attempt to estimate the strength of 
the magnetic field in these earlier systems, given some in- 
dependent estimate of the column density of free electrons, 
A^e- Following Kronberg and Perry (1982), we can express 
an RRMc arising at redshift zc as 

RRMc = /5(1 + zc)"^^e(S||)rad m"^ (13) 
where is in cm^'^, (Bn) accommodates any field reversal 



pattern, and is defined by 



J neB\\dl 

cloud 

J riedl 

cloud 



(14) 



The constant (3 — 2.63 xlO~^^ rad m~^ cm^ Gauss^^. 

The rms deviation from zero of the observed RRM, 
o'RRM, can be straightforwardly converted to an estimate 
of a typical by inverting Eq. (fT3|) and inserting an 

estimate of Ne for the system generating the Faraday ro- 
tation. The important point is that as we move into the 
redshift range 1-8 < z < 3.7 the fraction of lines of sight, 
fi, with near zero RRM, falls off sharply. This is indepen- 
dent of whether the excess RRM is in the vicinity of the 
radio source ("intrinsic"), or located in a galaxy system 
at intervening redshift. Note that the exact redshift is not 
crucial for our estimate here as the spread in (1 -I- z)^ over 
the interval 1-8 < z < 3.7 is modest, of order a few at 
most. 

Inverting Eq. (fT3|l and setting RRMc ^ ctrrm m the 
observer's frame (e.g. panel 9 of Fig[3l) gives 



B«) = 5.5 X lO^^G 



l + zc 
3.5 



fRRM 



20 rad m 



1.7 X 102icm-2 



(15) 



for the system causing the RRM excess. 

The local < -B||> values within a system are expected 
to be larger because of cancelations due to reversals in the 
magnetic field orientation along the line of sight. The col- 
umn density assumed in the above normalization is com- 
parable to HII column densities through today's typical 
spiral galaxies. Thus if the RRM is due to galactic systems 
at high redshifts, our estimate would imply the presence of 
magnetic fields there with strength of order at least a /xG, 
which is comparable to the value observed in the Galaxy. 
If the radio source is embedded within a Faraday rotating 
cloud (the intrinsic RRM case), the ~ 2 x lower RRM 
pathlength would, relative to a lower z intervening cloud, 
make the magnetic field strength correspondingly higher. 
The results presented in §4.1 indicate that magnetic fields 
must have been generated very quickly in galaxy systems, 
at early cosmological epochs. 
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5 10 15 , 20 25 5 10 15 , 20 25 

T.J2 (rad m-2) T.J2 (rad m"^) 

Fig. 8. — Contours of constant likelihood as a function of fit parameters V jlintr and cTdoud for two fixed values of <T„oise- Mgll systems 
with Wr > O.SA are used as a model for the intervening systems. 



6. ON THE DETECTABILITY OF A WIDESPREAD IGM 
MAGNETIC FIELD 

Limits on a widespread IGM magnetic field were first 
discussed and derived in the 1970's under the assumption 
that the Universe was homogeneous, had only baryonic 
matter, and was pervaded by contiguous magnetized cells 
of comoving length, Iq. Under these assumptions, RM's 
measured at the largest redshifts (then ^ 2) placed upper 
hmit estimates of B^'^^ in the range 10~® - 10"^ G (Rees 
& Reinhardt (1972), Nelson (1973), Kronberg & Simard- 
Normandin (1976), Kronberg et al. (1977)). In this co- 
expanding model in which |i?(z)| oc (1 + z)^, the observed 
growth in RRM(z) due to a widespread, co-expanding 
IGM is dominated by contributions at the largest redshifts. 
In the context of current concordance cosmology, and our 
knowledge of large scale filaments and voids of matter, 
these Bl'^^ limits are no longer appropriate, and need to 
be revisited. 

Now, isolating or limiting a widespread extragalactic 
RRM, and hence B^'-^^, is more difficult. It requires 
that we disentangle a B^'^^ from any cosmic evolution in 
RRMc(z) (case 2) and RRMintr (case 3). At least some of 
the high redshift sources in our sample, such as 3C191 at 
z=1.95 (Kronberg et al. 1990), and several high-z RM's 
measured by CariUi et al.(1997), Athreya et al. (1998), 
and Pentericci et al. (2000) indicate evolutionary effects 
in RRM„t,.. 

However the expected width of N(RRMigm, z) due 
to a widespread, co-expanding magnetized IGM increases 
steeply at the large redshifts (e.g. Kronberg, Reinhardt & 
Simard-Normandin 1977). This raises the possibility that 
future RRM datasets that extend to redshifts z = 4 - 5 
may be strongly influenced by a widespread co-expanding 
intergalactic medium. 

Then with better knowledge of the intervenor and 
source-intrinsic populations, RRMc(z) and RRMi„ir(^), 
it may be possible to explore, or limit the contribu- 
tion of RRMigm(z)- The data could then be compared 
with modelled contributions of RRMigm(2) (filaments) 
and RRMiGM(2)(voids) based on LSS evolution simula- 
tions. 

A probe of the RRM contribution of some local uni- 
verse galaxy filaments was recently attempted by Xu et 
al. (2006). Their initial estimate was -BjGM ^ 3 X 
10^''G in the Perseus-Pisces supercluster filament zones, 
scaled to an assumed field reversal scale of ^ 0.5 Mpc 
and to estimates of the electron density in the warm-hot 
IGM(WHIM). 



7. SUMMARY AND CONCLUSIONS 

We have analyzed an unprecedented large sample of 
RRM data extending to redshift 3.7. The size of our new 
Faraday RM sample permits us to isolate preferred Galac- 
tic {l,b) zones where Galactic foreground RM is optimally 
removed with the help of a new, more accurate set of 1566 
extragalactic source RM's that are mostly off the Galactic 
plane. We find a clear increase with redshift of Galaxy- 
corrected RRM's at medium to low RRM levels, below 
~ 100 rad m^^ in the observer's frame. Our results and 
conclusions can be summarized as follows: 

1. There is a striking and systematic decline, signifi- 
cant at the ~ 3 cr level, in the fraction of sources in 
the range z ^ 1.5 - 2.3 that have |RRM| less than 20 
rad m~^. This is the first global indication that the 
Universe becomes increasingly "Faraday opaque" to 
the highest redshift radio sources. 

2. A model with intervening systems distributed ac- 
cording to the dN/dz statistics of strong Mgll ab- 
sorption line systems in QSO spectra is consistent 
with the observed growth in the widths of the RRM 
distribution up to z 2 with no evolution in the 
rest frame RM of each absorber. This hypothesis 
could be tested in the future by comparing the dis- 
tribution of RM in quasar sightlines that have, and 
do not have, strong Mgll absorption. 

3. We provide global estimates of early Universe mag- 
netic field strengths in galaxy systems, based on 
the progressively increasing "Faraday opaqueness" 
of systems beyond z ~ 1.5. These have at least 
fiG level fields, and indicate that magnetic fields at 
these high redshifts were at least as strong as at the 
present epoch. 

4. Redshift variations in the widths of RRM distri- 
butions below z ~ 1.2 are small and can only be 
specified at the 1.5 - 2a level, limited by residual 
uncertainties in the foreground Galactic RM. The 
latter is at the level of 9 rad m^^ in the higher b 
zones selected, determined from a new all-sky sam- 
ple of 1566 RM's. 
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